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I 



We investigate an optically driven quantum computer based on electric 

> 

■ dipole transitions within coupled single-electron quantum dots. Our quan- 
tum register consists of a freestanding n-type pillar containing a series of 

o 

■ pair wise coupled asymmetric quantum dots, each with a slightly different 
"j^ ■ energy structure, and with grounding leads at the top and bottom of the pil 

lar. Asymmetric quantum wells confine electrons along the pillar axis and 
a negatively biased gate wrapped around the center of the pillar allows for 
electrostatic confinement in the radial direction. We self-consistently solve 



o 



^ . coupled Schrodinger and Poisson equations and develop a design for a three- 

qubit quantum register. Our results indicate that a single gate electrode can 
be used to localize a single electron in each of the quantum dots. Adjacent 
dots are strongly coupled by electric dipole-dipole interactions arising from the 
dot asymmetry, thus enabling rapid computation rates. The dots are tailored 
to minimize dephasing due to spontaneous emission and phonon scattering 
and to maximize the number of computation cycles. The design is scalable to 
a large number of qubits. 

PACS Number (s): 03.67.Lx, 73.20.Dx, 85.30.Vw 
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I. INTRODUCTION 



The possibility that a computer with exceptional properties could be built employing the 
laws of quantum physics has stimulated considerable interest in searching for useful algo- 



and factorization,a have been discovered; others, including the suggestion that quantum 
computers will prove useful to model quantum systems, are being sought. Meanwhile, var- 
ious physical implementations are being explored, including trapped ions,li cavity quantum 
electrodynamics,0 ensemble nuclear magnetic resonance,! small Josephson junctions,! optical 




based on quantum dots. Although the advantages of quantum computing are enormous 
for particular key applications, the requirements for their implementation are extremely 
stringent, perhaps especially rigorous for solid-state systems. Nevertheless solid-state quan- 
tum computers are very appealing relative to other possible implementation schemes because 
of the well-known ability to customize the design through the use of artificially structured 
materials and the probable scalability of the resulting design. For example, integrated circuit 
manufacturing technology would be immediately applicable to quantum computers of the 
proper implementation; and such designs would not only be scalable to smaller dimensions 
along the "semiconductor learning curve" but also large ensembles of "identical" quantum 
computers could be manufactured, that could be individually fine-tuned electrically. To 
date, no solid-state implementation of quantum computing has been demonstrated. 

In this paper, we investigate a solid-state quantum computer implementation that 
is amenable to manufacturing with integrated circuit technology. We develop a three- 
dimensional (3D) device model and self-consistently solve coupled Schrodinger and Poisson 
equations to generate a quantum computer design for a three-qubit quantum register that is 
based on pair wise coupled asymmetric III-V quantum dots. The design is optimized for a 
long coherence time and a rapid computation rate. Our results indicate that this structure 
may provide a realistic scalable candidate for quantum computing in solid-state systems. 
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II. PROPOSED STRUCTURE 



The proposed quantum dot quantum computer (see Fig. |^) consists of a pillar structure 
composed of a chain of asymmetric quantum dots separated by intervening layers of higher 
bandgap composition fabricated in a GaAs/AlGaAs technology by means of a sequence of 
planar MBE growth steps and subsequent etching to form the pillar. A sheath of similar 
AlGaAs composition is then grown surrounding the pillar and a wrap-around gate electrode 
deposited. A drain (source) is formed at the top (bottom), the series of asymmetric quan- 
tum dots are in the center region, and the gate surrounds the region of the pillar containing 
the quantum dots. Tarucha et al.El have reported similar n-type single electron transistor 
(SET) structures. Electron confinement along the pillar axis is produced by the band gap 
discontinuity of the dot structure. Encasing the quantum dot structure in the pillar core by 
the cylindrical sheath and the gate electrode provides confinement in the radial direction. 
By applying a negative bias that depletes carriers near the surface, an additional parabolic 
electrostatic potential is formed that allows for tuning of the radial confinement and local- 
ization of one electron per dot. The simultaneous insertion of a single electron per dot is 
accomplished by lining up the quantum dot ground state levels so that they lie close to the 
Fermi level; a single electron is confined in each dot over a finite range of the gate voltage 
due to shell filling effects.^ 

Thus, the pillar consists of a vertical stack of coupled asymmetric GaAs/AlGaAs quan- 
tum dots of differing size and composition so that each dot possesses a distinct energy 
structure. Qubit registers, |0) and are based on the ground and first excited state of 
the single electron within each quantum dot. Overall, parameters of the structure can be 
chosen to produce a well-resolved spectrum of distinguishable qubits. The asymmetric dots 
produce large built-in electrostatic dipole moments between the ground and first excited 
state, and electrons in adjacent dots are coupled through the electric dipole-dipole inter- 
action, while coupling between non- adjacent dots is significantly weaker. This produces 
the desired quantum computer consisting of a linear array of binary states (qubits) with 



3 



pair wise pillar-axis coupling between adjacent qubits.0 In addition to energy tuning, the 
asymmetry of each quantum dot can be designed so that dephasing due to electron-phonon 
scattering and spontaneous emission is minimized. The combination of strong dipole-dipole 
coupling and long dephasing times make it possible to perform many computational steps 
before loss of coherence, in fact, it is believed possible to design this device so that error 
correction substantially prohibits coherence loss. 

Quantum computations are performed by means of a series of coherent optical pulses in 
the far infrared. Final readout of the amplitude and phase of the qubit states can be achieved 
through quantum state holography. Amplitude and phase information are extracted through 
mixing the final state with a reference state generated in the same system by an additional 
delayed laser pulse and detecting the total time- and frequency-integrated fluorescence as a 
function of the delay.lli Extracting the final state information using quantum state holog- 



raphy requires multiple experiments, one for each delay, as described in Ref. |T^. Thus, the 
computation must be performed several times before an answer is arrived at. This is no 
real problem since the number of repetitions needed is only on the order of 40 or so, inde- 
pendent of the number of computational steps in a given quantum algorithm. Through the 
use of integrated circuit manufacturing technology, it is possible to simultaneously fabricate 
a large array of "identical" pillar quantum dot quantum computers, that is, on the order 
of 10^'' per wafer. Each of these quantum registers could be electrically connected through 
deposited interconnect in such a manner so that each could be individually tunable to pro- 
duce an array of identical units. In general, inhomogeneity among the quantum dots will 
result in slightly different energy levels. Sherwin et al.0 have recently pointed out that one 
can perform accurate qubit operations in an inhomogeneous population of quantum dots 
arising from quenched disorder due to static charged defects, for example, provided that 
each SET is independently calibrated. This calibration can done by performing simple gate 
operations and tuning the gate electrodes appropriately. Efficient optical coupling to the 
resulting ensemble can be achieved through optical light guiding as suggested in Ref. |T7[ By 
this means direct observation of fluorescence is possible. Quantum computations are per- 



formed by means of a series of coherent optical pulses in the far infrared, and may be carried 
out in complete analogy with the operation of an NMR quantum computer.0 It should be 
remarked that while our scheme resembles NMR ensemble quantum computation in the use 
of a series of optical pulses to perform quantum logic gates, it differs from NMR quantum 
computation in that our use of a collection of single electron transistors is done to enable a 
stronger signal to noise ratio in the readout phase. In principle, the quantum computation 
could be done with only a single SET transistor structure if the readout measurements were 
sufficiently sensitive. 



III. DEVICE MODEL 

In the context of studies of the Coulomb blockade in self-organized quantum dots and 
planar single-electron transistors, self-consistent calculations of electronic structure, shell 
filling effects, electron-electron interaction. Coulomb degeneracy, and Coulomb oscillation 
amplitudes have been carried out for various quantum dot structures.&S Our quantum 
register can be analyzed using methods similar to those used to study the self-consistent 
electronic structure in single-electron transistors. The problem we address is similar to those 
addressed by other authors who are interested in obtaining current-voltage characteristics 
and studying Coulomb oscillations in single-electron transistors over a wide range of gate 
biasing and shell filling conditions.EHHl 

In our case, we are interested in obtaining the self-consistent electrostatic potential and 
electronic eigenstates in an equilibrium configuration in which the source and drain are 
grounded and the gate electrode is negatively biased. The electrostatic potential, V{r), is 
obtained by solving the Poisson equation for n-doped semiconductorsil'0 

4-TT r n 

VV(r) = -— g [n{r) + N+{r)] (1) 

In the Poisson equation, q is the absolute value of the electron charge, e is the static dielectric 
constant, n(f^ is the electron concentration and N^{r) is the known concentration of ionized 
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donors in the structure. For the dielectric constant, we adopt the GaAs value e = 12xn The 
Poisson equation is solved subject to boundary conditions on the electrostatic potential, 
V{r). At the interfaces between the semiconductor and the source, drain and gate electrodes, 
V{r) is equal to the applied gate voltage while at the semiconductor- vacuum interfaces, the 
normal derivative of V{r) vanishes. 

Following Ref. the global electron concentration, n(r), in the device is obtained by 
partitioning the pillar structure into "bulk" and "quantum" regions. In the "bulk" regions 
far from the quantum wells i.e. the source and drain regions, electrons are treated in the 
Thomas-Fermi approximation and the electron concentration is given byil 



1 
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n{r) = \ ^ (2) 
otherwise 



where /i is the chemical potential and Vif) is the effective electron potential. The chem- 
ical potential, /x, is determined through the requirement that overall charge neutrality be 
maintained in the bulk regions, i.e. the chemical potential is adjusted until 

n{f)-N^{f) ) df=0 (3) 

where the integration is carried out over the bulk source and drain regions. 

The effective potential, f/(r), in the bulk regions includes the Hartree potential, Uh = 
—q V{r), and the conduction band offset, AE'c, which depends on the local Al concentration, 
X. Thus, 

U{f) = -q V{r) + AE^ (4) 

where the conduction band offset, AEc, is taken to be 60% of the difference between the 
AlxGai-xAs and GaAs bandgaps. Using the bandgap variation of Al^Gai-xAs determined 
by Lee et al.,0 we obtain the following expression for the conduction band offset as a function 
of the local Al concentration, x: 

AE^ = 0.6 (1155 X + 370 x^) meV (5) 



In the "quantum" regions containing the quantum wells, the electron concentration, n(r), 
is determined by the electron wavefunctions, ipi{r), and energies, Ei, through the relation 

'^(^ = H'^i (6) 
i 

The electron occupancy in each level, n^, is a function of the electron energy and the tem- 
perature. 

The electron wavefunctions and energy levels, E^, are obtained by solving the Schrodinger 
equation in the effective mass approximation 

^i(r) = (7) 

The electron potential, t/(r), in the quantum regions is given by 

U{f) = -q V{r) + AE, + U,,{r) (8) 

where Uxci,'^) is the exchange-correlation potential of Perdew and Zunger.0 

In the quantum register discussed in the next section, the gate voltage is negatively 
biased in such a way that a single electron is strongly localized in each electrostatically 
confined quantum dot. The radial confinement potential is strong enough that the lowest 
few electron wavefunctions are strongly localized near the center of the pillar and die away 
far from the semiconductor-electrode interface. In our design, the quantum wells are wide 
enough and the barriers between the quantum wells are thick enough so that the lowest few 
electron wavefunctions do not penetrate to the center of the barriers separating the quantum 
wells. Since all the wavefunctions of interest vanish at the center of these barriers, we can 
divide the quantum well region into several regions (one for each qubit). These regions are 
taken to be cylinders stacked along the pillar axis with top and bottom surfaces located at 
the centers of the barriers between adjacent wells. We solve the Schrodinger equation in 
each dot separtely subject to the boundary condition that all wavefunctions vanish at the 
region boundaries. 

Due to the cylindrical symmetry of the structure, the 3D Schrodinger equation can 
be reduced to a 2D equation in cylindrical coordinates. One might try to solve the 2D 
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JL^V^ + U{t) - E, 
2 m* 



Schrodinger equation by finite differencing the partial differential equation and solving the 
resulting matrix eigenvalue equation. The size of the matrix to be diagonalized is equal to the 
number of interior mesh points in the 2D grid and this is much too large to be handled easily. 
Other authors have taken this brute-force approach to solving the Schrodinger equation in 
self-consistent Poisson-Schrodinger problems with the result that solving the Schrodinger 
equation is the most time consuming part of the computation.ii We find that it is possible 
to do better. In solving the 2D Schrodinger equation, we first approximate U{p, z) in each 
quantum dot by a separable potential 

U{p,z)^Us{p,z) = Ur{p) + U,{z) (9) 

where the axial potential is defined as 

U.iz) = ^j\{p,z)pdp (10) 

and the radial potential is given by 

Urip) = ]; l\u{p,z)-U,{z) ) dz (11) 
L Jo 

In these last two expressions, R and C are the radius and height of the cylindrical region 
over which U{p, z) is defined in each dot. With the separable potential approximation, 
the 2D Schrodinger equation can be separated into two ID equations which can be cast as 
finite difference eigenvalue equations and solved numerically for the electron energies and 
wavef unctions. The resulting 2D wavefunctions are the best product wavefunctions that ap- 
proximate the solution of the 2D Schrodinger equation in each qubit. The electronic states 
in the separable potential approximation in our cylindrical pillar are labeled by three quan- 
tum numbers (np, n^, n^) which specify the number of nodes in the product wavefunctions 
associated with cylindrical coordinates p, (p, and z. In this notation, the qubit state |0) is 
denoted (0,0,0) while |1) is denoted (0,0, 1). We find that the separable wavefunctions are 
reasonable approximations to the true wavefunctions since we are starting with a separable 
potential which is already close to the true potential in some average sense. We next obtain 
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the exact energies and wavefunctions of the original non-separable Schrodinger equation by 
treating the residual U{p, z) — Us{p, z) as a perturbation and expanding the exact wavefunc- 
tions as a sum of separable wavefunctions. Our expansion of the true wavefunctions in terms 
of separable wavefunctions is rapidly converging and we find that the dominant terms in the 
expansion of the true wavefunctions are the separable wavefunctions of the same symmetry. 
Our approach to solving the 2D Schrodinger equation is fast and most of computing time is 
spent solving the Poisson equation. 

To complete the specification of the electron charge density in the quantum dots, it is 
necessary to compute the electron occupation numbers, rij. One might expect that rii would 
be given by the Fermi-Dirac distribution and indeed this would be the case if the electrons 
in the dots were delocalized and in tunneling contact with the leads. In this case, the qubits 
could exchange electrons with their environment and the total number of electrons in the dot 
N = J^ifii could take on non-integer values. But clearly this is not tolerable in a quantum 
computer and we must carefully arrange things so the dot wavefunctions exhibit a high 
degree of localization. In this situation, only an integer number of electrons can occupy the 
dot and this constraint gives rise to what is known as the Gibbs distribution. The number 
of electrons, N, in the dot is determined by minimizing the Gibbs free energy with respect 
to the integer number of electrons, N. The Gibbs free energy is F{N) = —kT \n[Z{N)], 
where the grand canonical partition function, Z{N), is given byi'i 



The lack of diffusive contact between the quantum dots and the rest of the device means 
that the chemical potential, /z, is determined by electrons in the leads and contacts. The 
summation in Z{N) is carried out over all electron configurations {rii} for which J2i = N. 
Double counting the Coulombic interaction is avoided by subtracting the Hartree energy 
Eh{N) for the N electrons. The Hartree energy appearing in the partition function is@ 



Z{N) = Y: exp 

{rii} 



Y.^niE,-EH{N)-^lN 
kT 



(12) 




(13) 
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where ne{r) is the charge in the quantum dot and the integration is restricted to the dot 
region. Directly solving for the Hartree energy by performing a double integral over the 
quantum dot charge density is too time consuming and impractical due to the presence of 
the singularity in the integrand at r*= r*. An alternative method of calculating the Hartree 
energy is to use the equivalent expression 



where the potential Ve{r) is obtained by solving the Poisson equation in the pillar using the 
charge density, ne{r), in the quantum dot^ 

The boundary condition on Ve(r) at the surface of the pillar is determined by asymp- 
totically expanding Ve{r) in a multipole expansion in the quantum dot charge density up 
through quadrupole terms and using this expansion to specify Ve at the surface. This is 
a good approximation since the pillar boundaries are far from the localized quantum dot 
charge.0 

To obtain a self-consistent solution to the coupled Poisson and Schrodinger equations, 
we first specify the device structure including the Al^Gai-xAs alloy composition, the doping 
concentration, and the arrangement of the electrodes. In all our runs, the source and drain 
are assumed to be grounded and the gate is assumed to be negatively biased. We initially 
assume complete depletion in the structure and solve the Poisson equation to obtain an 
initial guess for the electrostatic potential, V{p, z). With this electrostatic potential and the 
quantum well band offset potentials, we solve the Schrodinger equation for the unoccupied 
quantum dot energies and wavefunctions. The chemical potential in the depleted structure 
is set to the minimum of the Thomas- Fermi electron potential, U = —q V{p, z)+AEc, in the 
source and drain regions. Starting with these initial guesses for the chemical potential in the 
leads and the solutions of the Poisson and Schrodinger equations, we obtain self-consistent 
solutions through the following relaxation procedure. 

First, electron densities in the leads and the quantum dots are obtained from the chemical 
potential, the temperature, and the quantum dot electronic states. The global charge density. 




(14) 



10 



including the given doping charge, is then obtained. Next the Poisson equation is solved 
for V{p,z). The Hartree potential and exchange-correlation potentials are then obtained 
from V{p, z) and the electron charge density. With the electron potentials in hand, the 
Schroodinger equation is solved in each quantum dot region. The procedure is then repeated 
until convergence is achieved. In updating the electrostatic potential and electron charge 
density, the new solutions are mixed with the old to obtain the updated solutions. For the 
electrostatic potential 

V{p, z)^X V;e^(p, z) + {l- A) KwIp, z) (15) 

where A < 1 is a relaxation parameter which is dynamically adjusted to accelerate conver- 
gence. A similar scheme is used to update the electron charge density. The above procedure 
is iterated until the chemical potential, electrostatic potential, electron charge density, and 
quantum dot energy levels all change by less than some small relative tolerance between 
successive iterations at which point convergence is achieved. Typically about 400 iterations 
are required to achieve convergence to within one part in 10^. 

IV. A THREE QUBIT QUANTUM REGISTER: ID ANALYSIS 

We can use the device modeling program described in the last section to obtain a design 
for a three-qubit quantum register. We could, in principle, do a full 3D analysis of the 
device and obtain suitable design parameters (i.e., pillar dimensions, doping concentrations, 
asymmetric well shapes, electrode placement and biasing, etc.) based on our computationally 
intensive 3D model. Clearly this would be prohibitively time consuming due to the size of 
the parameter space that would need to be investigated as well as the time required to 
perform each run. To narrow down the design parameters, we can take advantage of the 
fact that our quantum computer is operated in the extreme depletion limit and do a simple 
ID analysis to gain some useful insight. 

Let's assume that inside the core of stacked quantum wells (radius Rc) we have complete 
depletion and uniform doping. In this limit, the quantum dot electron potential, U (r), can be 
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expressed in cylindrical coordinates as U{r) = U{z) + U{p), where U{p) is a radial potential 
arising from the uniform donor density and U{z) is the conduction band offset potential 
along the growth direction. This separable potential assumption is a good approximation in 
the strong depletion regime where only a single electron resides in each dot. The assumption 
of a separable potential is commonly used in the study of quantum dot structures and enables 
us to consider the z and p motions separately.El'ii The z-directional potential U{z), shown 
schematically in the inset of Fig. ^, is a step potential formed by a layer of Al^Gai-xAs 
of thickness B {0 < z < B) and a layer of GaAs of thickness L — B {B < z < L). The 
resulting asymmetric quantum dot/well is confined by AlyGai^yAs barriers with y > x and 
the asymmetry is parameterized by the ratio B/L where < B/ L < 1. 

In the effective mass approximation, the qubit wavefunctions are \i) = R{p) ipiiz) Us{r) 
{i = 0, 1). Here R{p) is the ground state of the radial envelope function, ilJi{z) is the envelope 
function along z, and Us{r) is the s-like zone center Bloch function including electron spin. 
For simplicity, we assume complete confinement by the AlyGai^yAs barriers along the z 
direction. Then, the envelope function ipi{z) is obtained by solving the time-independent 
Schrodinger equation subject to the boundary conditions = 'ipiiL) = 0. The energies 

of the qubit wavefunctions are given hj E = Ep + Ei where Ep is the energy associated with 
R{p) and Ei is the energy associated with ilJi{z). 

Figure H shows the probability density, as a function of position, z, for the two 

qubit states |0) and |1) in a 20 nm GaAs / Alo^^GaQ^-jAs asymmetric quantum dot. The barrier 
thickness 5 = 15 nm and the overall length of the dot is L = 20 nm. By choosing B/L = 
0.75 and x = 0.3, it is found that the ground state wavefunction |0) is strongly localized in 
the GaAs region while the |1) wavefunction is strongly localized in the Alo.3Gao.7As barrier. 
By appropriately choosing the asymmetric quantum dot parameters, the qubit wavefunctions 
can be spatially separated and a large difference in the electrostatic dipole moments can be 
achieved. 

The transition energy AEq = Ei — Eo between |1) and |0) is shown in Fig. ^ as a 
function of B/L in a 20 nm GaAs/Alr^Gai-xAs asymmetric quantum dot {L = 20 nm). 

12 



Several values of Al concentration x are considered. In Fig. ^ we fix the Al concentration at 
X = 0.2 and plot AEq as a function of L for several values of B/L. The continuous curves 
are based on our ID analysis and the squares are the qubit energy gaps for a three qubit 
self-consistent quantum register calculation as described in the next section. It is clear that 
the transition energy can be tailored substantially by varying the asymmetry parameter. 
With three parameters available for adjustment {B, L, and x), we can make AEq unique for 
each dot in the register. In this way, we can address a given dot by using laser light with 
the correct photon energy. 

It is desirable that the |1) state be the first excited level of the quantum dot. Thus, 
the lowest lying radial state (0, 1,0) should lie above the |1) state. The radial energy gap, 
AEi, between the ground state, |0), and the first radial excited state, (0, 1,0), is found by 
solving a 2D Schrodinger equation for an electron in the radial potential, Ur{p). If we take 
the barrier in the sheath to be infinite, then in the extreme depletion limit, we have 



where V{p) is the radial electrostatic potential. For complete depletion and uniform doping, 
the Poisson equation for V{p) can be solved analytically. Thus, 



where Rc is the sheath radius and N]^ is the doping density. Numerically solving the 2D 
Schrodinger equation for an electron in the potential, Ur{p), is straightforward. Figure |^ 
shows the radial energy gap, AEi, between the |0) and the lowest lying radial state, (0, 1,0), 
as a function of doping concentration, N^, for several values of Rc- For narrow pillars 
with low doping concentrations, the radial energy gap is determined by size confinement. 
For large pillars with high doping concentrations the radial energy gap is determined by 
electrostatic confinement. From Fig. ^, we see that the qubit energy gaps reach a minimum 
near AEq ~ 70 meV for quantum wells with L ~ 20 nm. The results of Fig. ^ suggest that 




(16) 



V{p) = ^{Rl - p') 



(17) 
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radial gaps in the range of AEi ~ 100 meV with strong size confinement can be achieved 
with doping densities in the range of ~ 10^^ cm^^ if Rc ~ 70 A. 

The electric field from an electron in one dot shifts the energy levels of electrons in adja- 
cent dots through electrostatic dipole-dipole coupling. By appropriate choice of coordinate 
systems, the dipole moments associated with |0) and |1) equal in magnitude but oppositely 
directed. The dipole-dipole coupling energy is then defined sM 



\di\ \d' 



Vm = 2 (18) 



where di and d2 are the ground state dipole moments in the two dots, = 12.9 is the 
dielectric constant for GaAs, and Ru is the distance between the dots. 

Figure | shows the dipole-dipole coupling energy, Vdd, between two asymmetric 
GaAs/Al^Gai-x^s quantum dots of widths LI = 19 nm and L2 = 21 nm separated by 
a 10 nm AlyGai-yAs barrier. The coupling energy is plotted as a function oi B/L for sev- 
eral values of x where B/L and x are taken to be the same in both dots. The dipole-dipole 
coupling energies are a strongly peaked function of the asymmetry parameter, B/L. From 
the figure, we see that values of Vdd ~ 0.15 meV can be achieved. 

Quantum dot electrons can interact with the environment through the phonon field, par- 
ticularly the longitudinal-optical (LO) and acoustic (LA) phonons. The LO phonon energy, 
Tiujloi lies in a narrow band around 36.2 meV . As long as the quantum dot energy level 
spacings lie outside this band, LO phonon scattering is strongly suppressed by the phonon 
bottleneck effect. Acoustic phonon energies are much smaller than the energy difference, 
Ai?, between qubit states. Thus, acoustic phonon scattering requires multiple emission pro- 
cesses which are also very slow. Theoretical studies on phonon bottleneck effects in GaAs 
quantum dots indicate that LO and LA phonon scattering rates including multiple phonon 
processes could be slower than the spontaneous emission rate provided that the quantum dot 
energy level spacing is greater than ~ 1 meV and, at the same time, avoids a narrow window 
around the LO phonon energy^^ In Ref. ^ Inoshita and Sakaki compute multi-phonon 
relaxation rates in spherical single-electron GaAs quantum dots due to one- and two-phonon 
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scattering by LO and LA phonons at T = and T = 300 K. Using the results of this 
calculation, we estimate that multi-phonon scattering dominates the spontaneous emission 
only if the qubit energy level spacing is within ~ 4 meV of the LO phonon energy. Likewise, 
multi-phonon LA scattering becomes important if the qubit energy gaps are smaller than 
~ 1 meV 

While dephasing via interactions with the phonon field can be strongly suppressed by 
proper designing of the structure, quantum dot electrons are still coupled to the environment 
through spontaneous emission and this is the dominant dephasing mechanism. Decoherence 
resulting from spontaneous emission ultimately limits the total time available for a quantum 
computation.!! Thus, it is important that the spontaneous emission lifetime be large. The 
excited state lifetime, T^, against spontaneous emission iS 

where D = (0|2;|1) is the dipole matrix element between |0) and 

Figure |^ shows the spontaneous emission lifetime of an electron in qubit state |1) for a 
20 nm GaAs/Al^Gai^x^s quantum dot as a function of asymmetry parameter, B/L, for 
several values of Al concentration, x. It is immediately obvious from Fig. |^ that the lifetime 
depends strongly on B/L. Depending on the value of x chosen, the computed lifetime can 
achieve a maximum of between 4000 ns and 6000 ns. In general, the maximum lifetime 
increases with x. In Eq. (0), the lifetime is inversely proportional to AE^ and D^, but the 
sharp peak seen in Fig. ^ is due primarily to a pronounced minimum in D. 

Based on these results, we can estimate parameters for a solid state quantum register 
containing a stack of several asymmetric GaAs/AlQ^GaQ jAs quantum dots in the L ~ 20 nm 
range separated by 10 nm AlyGai^yAs barriers {y > 0.4). An important design goal is 
obtaining a large spontaneous emission lifetime and a large dipole-dipole coupling energy. 
From Figs. ^ and |^, we see that both can be achieved by selecting an asjTumetry parameter, 
B/L = 0.8. This gives us a spontaneous emission lifetime Td = 3100 ns and a dipole-dipole 
coupling energy Vm = 0.14 meV. The transition energy between the qubit states is on 
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the order of 100 meV (A = 12.4 fim). In a quantum computation, the quantum register 
is optically driven by a laser as described in Ref. In our example, we require a tunable 
infra-red laser in the mid- 10 fim range so we can individually address various transitions 
between coupled qubit states. 

V. A THREE QUBIT QUANTUM REGISTER: 3D ANALYSIS 

Using the results of our simple ID model as a starting point, we designed a three qubit 
quantum register by using the self-consistent device model described in Section III. Several 
criteria have to be met for a viable quantum register design and the structure we obtained 
through trial-and-error involved tradeoffs between several design goals. 

For a self-consistent quantum register calculation, we assume the parameters of the free- 
standing quantum dot pillar structure (shown in Fig. |1|) as follows: The height of the pillar 
is taken to be £ = 1000 nm while the radii of the core and sheath are taken to he Rc = 7 nm 
and R = 50 nm. The drain and source contacts at the top and bottom of the pillar are 
grounded and a cylindrical gate with a height of 400 nm is placed around the center of the 
pillar. Near the source and drain contacts, layers of intrinsic semiconductor serve to inhibit 
gate-to-source and gate-to-drain currents. The central 600 nm of the pillar is uniformly 
n-doped with a doping concentration of A^^^ = 5 x 10^^ cm~^ . 

The cylindrical sheath surrounding the core region is composed of high band gap 
y4/o.45Gao. 55^-5 and serves to confine electrons to the core region. The three qubits in the 
core are defined by the composition profile of Al^Gai^x^s along the pillar axis. In our 
structure, the Al concentration, x, in the core region is uniform in the radial direction. The 
composition profile along the pillar axis in the core region is shown in Fig. The ground 
and first excited electronic states are the qubit states |0) and |1) and the electron charge 
densities for these states are shown schematically in the figure. We find that in thermal 
equilibrium the electrons reside entirely in the ground state |0) for temperatures as high as 
77 K since the energy gap between |0) and |1) is much greater than kT. This is indicated 
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schematically by the solid circles in the diagram. The widths L of the asymmetric quantum 
wells/dots defining qubits 1 through 3 are 19.0, 20.5 and 22.0 nm while the B/L ratios are 
0.670, 0.683 and 0.675 respectively. Our ID analysis suggests that asymmetry parameters 
in this range will result in long spontaneous emission lifetimes and strong dipole-dipole cou- 
pling between neighboring qubits. The asymmetric quantum wells are composed of GaAs 
and AlQ,2Gao,sAs layers and the barriers between the asymmetric dots/wells are composed 
of AloA5Gao,55As. 

With a properly chosen reverse gate bias, Vg, the doping region in the center of the pillar 
is depleted and the equilibrium Fermi level lines up so that there is exactly one electron 
in each dot. Single electron occupancy in the dots is necessary in order for there to be 
a well defined quhit Hilhert space. Due to shell filling effects, single electron occupancy in 
all three dots holds over a finite range of the gate voltage. By running our device model 
for several values of V^, we find that single electron occupancy is obtained over the range 
— 1.56 V < Vg < — 1.48 V . Thus, the requirement for single electron occupancy in 
the quantum dots is maintained in the presence of gate voltage fluctuations on the order 
of AVg ~ 0.08 V . For Vg = —1.5 V, the self-consistent electron potential along the pillar 
axis, (i.e. p = 0) is shown in Fig. ^ as a function of position along the pillar axis. The 
position along the pillar axis is measured from the drain contact at ^ = nm to the source 
contact at z = 1000 nm. Figure |^ is centered on the active region of the register containing 
the three quantum dots and the origin of the energy scale is chosen to be the equilibrium 
Fermi level. The total electron potential is approximately the sum of the self-consistent 
electrostatic Hartree potential and the conduction band offset potential, the self-consistent 
exchange-correlation potential being negligible. 

The self-consistent electron levels are obtained by solving the Schrodinger equation in 
the self-consistent potential shown in Fig. ^. In our structure, the |0) ground states have 
(np, ?T,0, Uz) = (0, 0, 0) symmetry and the |1) states (the first excited level) in all three qubits 
are (n^, n,^, n^) = (0,0,1) states. The self-consistent qubit energy gap, AEq, between the 
|0) and |1) states, the radial energy gap, AEi, and the spontaneous emission lifetime of 
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the |1) state, r^, and dipole moment, d, for the three qubits are hsted in Table |. From 
Table |, we see that the radial energy gaps are larger than the qubit energy gaps. Another 
thing to note is that the qubit energy gaps are large compared to kT at T = 77 K. Thus, 
in thermal equilibrium the electrons reside entirely in the |0) level at 77 K. This means 
that the initial state of the quantum register is characterized by a pure state density matrix 
Po = |0,0,0) (0,0, 0|. Consequently, there is no need for initial state preparation in our 
quantum register. In Fig. |lO|, the self-consistent electron probability densities in the three 
quantum dots are plotted as a function of position along the pillar axis. Each dot traps 
one electron and the probability densities in the ground and first excited states are shown 
as solid and dot-dashed lines, respectively. The barriers are thick enough so that electron 
wavef unctions in adjacent dots do not overlap. 

The energy levels for the three qubit quantum computer are shown in Table |I| with and 
without the inclusion of dipole-dipole coupling between the qubits. From Table |I| we see 
that a different energy corresponds to each three-electron state 12-, is) of the register 
where i„ = (0, 1) labels the state of the n-th qubit. Transition energies between the states 
|0) and |1) for a given qubit are obtained by taking differences between the appropriate 
entries in Table ^ For the first qubit, we take differences between all three-particle states 
|0, ^2, is) and |1, i2, is)- In general, the transition energy between |0) and |1) for an 
electron in the first qubit will depend on the states, ^2 and is, occupied by the second and 
third qubits, and there can be as many as four such conditional transitions. In the absence 
of dipole-dipole coupling between qubits, all four conditional transition energies between |0) 
and |1) for a given qubit are degenerate. When dipole-dipole interactions between the qubits 
are considered, the four-fold degenerate conditional transition energies split into multiplets 
depending on which states are occupied by the electrons in neighboring qubits. 

The conditional transition energies between |0) and |1) states for our three qubit register 



are shown in Fig. |Tl] as a function of photon energy. In the absence of dipole-dipole coupling, 
the transition energies for the three qubits are 40.86 meV, 47.14 meV, and 52.31 meV, 
respectively. When dipole-dipole interactions between qubits are taken into account, the 



conditional transition energies split into multiplets as shown in this figure. Each transition 
in the spectrum is labeled by the neighboring electron states which give rise to it. By 
performing optical vr-pulses at selected conditional transition frequencies, quantum logic 
operations can be performed. For example, a tt pulse performed on the lowest energy 
transition in Fig. |TT| performs a bit fiip on the first qubit provided the second qubit is in 
state This operation is just a ControUed-Not gate with qubit 2 as the control bit and 
qubit 1 as the target bit. 

The need to selectively perform vr-pulses at the conditional transition frequencies allows 
us to make some preliminary estimates on the parameters of the laser system needed to 
drive a quantum computation. If we want to selectively drive a given transition without 
exciting neighboring transitions, then the bandwidth of the vr-pulse needs to be less than 



the splitting between the two most closely spaced lines in the spectrum. From Fig. [TT], the 
two most closely spaced lines are spaced Ahuj ^ 0.0776 meV apart. If we require that 
the TT-pulse bandwidth is AE^^ ^ 0.01 meV, then the pulse length can be estimated from 
Heisenberg's uncertainty principle, AE'^AT^ /i/2, to be T^r ~ 33 ps. If we assume a square 
TT-pulse, the magnitude of the optical electric field is given byil 

= (20) 

q dT^ 

where d is the optical dipole from Table | and the average Poynting vector during the pulse 



Sav = -ceoEl (21) 
For d ^ 10 A and ^ 33 ps, we obtain Eq ^ 0.627 kV/cm and Sav ~ 522 W/cw?. 



VI. SUMMARY 

In this paper, we have studied a solid state implementation of quantum computing based 
on coupled quantum dots. Our quantum register consists of a free standing n-type pillar 
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with grounding leads at the top and bottom of the structure. Asymmetric quantum wells 
confine electrons along the pillar axis and a high bandgap AlGaAs sheath wrapped around 
the center of the pillar allows for confinement in the radial direction. The ground and first 
excited electronic states of the quantum dots act as qubit states |0) and |1), respectively. We 
have developed a 3D device model for a general SET structure containing several quantum 
dots. We self-consistently solve coupled Schrodinger and Poisson equations for the device 
and develop a design for a three qubit quantum register with asymmetric quantum dots 
tailored for long dephasing time and large dipole-dipole coupling between the dots. Our 
results indicate that a single gate electrode can be used to localize a single electron in 
each of the quantum dots. Adjacent dots are strongly coupled by electric dipole-dipole 
interactions arising from the dot asymmetry thus enabling rapid computation rates. 
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TABLES 

TABLE I. Self-consistent qubit energy gap, AEq, radial energy gap, AEi, spontaneous emission 

lifetime, r^, and optical dipolc moment, d, for a three qubit quantum register. 

Qubit No. AEo {meV) AEi (meV) Ts (ns) d (1) 



1 40.86 63.7 28000 11.7 

2 47.14 63.2 19000 11.5 

.3 52. .31 61.9 llOnO 11.4 



TABLE IL Register energies for a 


three qubit quantum ref 


;ister, with and without dipole-dipole 


coupling interaction between qubits. 






Register state 


Energy (meV) 


Energy (meV) 




uncoupled 


dipole-dipole coupled 


|0 0) 


-70.161 


-70.357 


|0 1) 


-17.852 


-17.833 


|0 1 0) 


-23.017 


-22.822 


|0 1 1) 


29.292 


29.273 


|1 0) 


-29.292 


-29.311 


|1 1) 


23.017 


23.213 


|1 1 0) 


17.852 


17.871 


ll 1 1) 


70.161 


69.966 
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FIGURES 

FIG. 1. Schematic illustration of the proposed quantum dot pillar structure. As an example, 
the structure is shown for a three qubit quantum register (i.e., with three quantum dots in the 
middle) . 

FIG. 2. Probability density along the confinement direction, z, for the qubit wavefunctions |0) 
(solid line) and |1) (dot-dashed line). The inset shows a schematic illustration of the conduction 
bandedge profile in the z direction. 

FIG. 3. Transition energy, A£^, between |0) and |1) in an L = 20 nm GaAs / AlxGai-xAs 
asymmetric quantum dot as a function oi BjL for several values of x. 

FIG. 4. The transition energy between the |0) and |1) qubit states as a function of the asym- 
metric quantum well width, L, for several values of the asymmetry parameter BjL and fixed 
barrier Al concentration x = 0.2. The continuous curves are based on a simple ID analysis and the 
squares are the qubit energy gaps for quantum dots as determined by the self-consistent calculation 
described in the text. 

FIG. 5. The energy gap between the |0) ground state and the lowest lying radial eigenstate 
(with {np,n^,nz) = (0,1,0) symmetry) as a function of the n-doping concentration, N^. The 
radial energy splitting for several pillar radii. Greater confinement of the radial wavefunction can 
be obtained by decreasing the pillar radius (size confinement) or increasing the doping concentration 
(electrostatic confinement). In the limit of large pillar radii, the radial energy splitting is determined 
by confinement in a parabolic electrostatic potential. 

FIG. 6. Dipolc-dipolc interaction between two asymmetric GaAs / AlxGai-xAs quantum dots 
of widths LI = 19 nm and L2 = 21 nm separated by an AlyGai^yAs barrier of width Wb = 10 nm. 
The coupling energy is plotted as a function oi B/L for several values of x. B/L and x are the 
same for both dots. 
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FIG. 7. Spontaneous emission lifetime for qubit state |1) in a GaAs / AlxGai-xAs quantum dot 
with L = 20 nm as a function oi B/L for several values of x. 

FIG. 8. Composition profile along the pillar axis for a three qubit quantum register. 

FIG. 9. Self consistent electron potentials as a function of position along the pillar axis in the 
active region of a three qubit quantum register. The total potential (solid line) is the sum of the 

Hartree potential (dot-dashed line) and the band offset potential. The fermi level (dotted line) is 
aligned so that exactly one electron resides in each quantum dot. 

FIG. 10. Electron density as a function of position along the pillar axis in the active region of 
a three qubit quantum register. The bias voltage is adjusted so that exactly one electron resides 
in the ground state of each quantum dot. In each dot, the probability density of the ground states 
(solid line) are shown along with the probability density of the first excited states (dot-dashed 
line). Depending on which state the electrons occupy, the electric dipole moment can point either 
left or right. 

FIG. 11. Conditional transition energies between the qubit states |0) and |1) as a function of 
photon energy for a three qubit quantum register. In the absence of dipole-dipole coupling, the 
transition energies are 40.86 meV, 4^1. \A meV, and 52.31 meV respectively. When dipole-dipole 
interactions between qubits are taken into account, the transition energies depend on the state of 
the adjacent qubits. 
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